Dielectron production in Ar+KCl collisions at 1.76A GeV 
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We present results on dielectron production in 40 Ar+KCl collisions at 1.7QA GeV. For the first 
time uj mesons could be reconstructed in a heavy-ion reaction at a bombarding energy which is 
well below the production threshold in free nucleon-nucleon collisions. The uj multiplicity has 
been extracted and compared to the yields of other particles, in particular of the (j> meson. At 
intermediate e + e~ invariant masses, we find a strong enhancement of the pair yield over a reference 
spectrum from elementary nucleon-nucleon reactions suggesting the onset of non-trivial effects of 
the nuclear medium. Transverse-mass spectra and angular distributions have been reconstructed in 
three invariant mass bins. In the former unexpectedly large slopes are found for high-mass pairs. 
The latter, in particular the helicity- angle distributions, are largely consistent with expectations for 
a pair cocktail dominated at intermediate masses by A Dalitz decays. 
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I. INTRODUCTION 

Lepton pairs emitted from the hot and dense colli- 
sion zone formed in heavy-ion reactions are appropriate 
probes for investigating in-medium properties of hadrons 
and, in general, the properties of hadronic bulk matter 
under extreme conditions. In fact, according to theory, 
there is even potential to detect new phases of nuclear 
matter in the laboratory by isolating their telltale sig- 
nals, among which are the direct decays of the short- 
lived vector mesons into e + e~ or p + u~ pairs Indeed, 
the electromagnetic current-current correlator, which en- 
ters into the virtual photon emission rate , depends on 
the properties of the strongly interacting medium and its 
constituents Q. 

Dilepton spectra measured by the CERES [4] and 
NA60 [|| experiments at CERN-SPS energies (40A - 
158 A GeV) pointed to a significant in-medium modifi- 
cation of the p meson spectral function, as signaled by a 
large additional yield (excess) of lepton pairs in the in- 
variant mass region below the p meson pole mass. At 
the much lower beam energies of 1 - 2A GeV the DLS @ 
collaboration at the Bevalac observed a similar dielec- 
tron excess over the "hadronic cocktail", i.e. the cock- 
tail resulting from meson decays in the late (freeze-out) 
phase of the collision. However, in contrast to the situa- 
tion at higher energies, for a long time this excess could 
not be satisfactorily explained by theoretical models and 
became the so-called "DLS puzzle". The baryon-rich 
regime probed at low beam energies obviously requires a 
different theoretical treatment than the pion-dominated 
SPS regime. 

The excess of electron pairs observed by DLS in the 
eta mass region has recently been re-investigated with 
the HADES 1 detector @] at SIS18 with carbon beams 
of 1 and 2 A GeV bombarding carbon targets @, 
These new data fully confirmed the earlier DLS result and 
thereby re-challenged theory to come up with a proper 
description of pair production at low energies. Compar- 
ing the excitation function of the excess pair multiplic- 
ity with the known 7r° and n meson production [10l - ll2l | 
revealed that the excess scales with bombarding energy 
very much like the pion yield does, but not at all like the 
eta yield. This finding already suggested baryonic reso- 
nances - and mainly the A(1232) - as possible source of 
the e + e~ excess yield. 

Besides the role played by baryon resonance decays, 
also a strong virtual bremsstrahlung contribution to the 
pair yield from mainly n+p interactions had been pre- 
dicted in various microscopic model calculations [1 3l — ll 6j j . 
Evidently, a good understanding of the pp — > ppe + e~ 
and np — > npe + e~ channels is required for a firm in- 
terpretation of dielectron emission from heavy-ion col- 
lisions. First data on dilepton production in nucleon- 
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nucleon collisions had again been obtained by DLS [l7j j . 
although with limited mass resolution and, at the lowest 
beam energies, limited statistics. Initial attempts to de- 
scribe those results were based on the soft photon approx- 
imation [l3j , followed later by more involved calculations 
using the one-boson exchange (OBE) approach (l4T4l6| . 
OBE models include contributions from a number of di- 
agrams which add up coherently, leading to subtle, but 
sizeable interference effects in the cross sections of both 
n+p and p+p reactions. Moreover, extending those re- 
sults to the description of heavy-ion reactions, e.g. in the 
framework of transport models, is a very difficult task 
which is at present not yet fully mastered. 

On the experimental side things moved again when 
HADES started to study p+p and d+p interactions at 
Ekin = 1-25 A GeV, i.e. just below the free w meson 
production threshold. The main goal of these experi- 
ments was to understand the n+p bremsstrahlung con- 
tribution to e + e~ production and to establish an experi- 
mental cocktail of dielectrons from " free" hadron decays 
at SIS18 energies. Using a deuterium beam, the "quasi- 
free" np — > npe + e~ reaction was therein tagged by the 
detection of a forward-going spectator proton. A com- 
parison of these data with our former 12 C+ 12 C result 
showed that pair production in the light C+C system 
can be understood as resulting from a mere superposi- 
tion of free N+N collisions 18]. Moreover, the excess 
pair yield in the C+C system, when normalized to the 
pion multiplicity, turned out to be largely independent 
of bombarding energy in the range of 1 to 2 A GeV, thus 
providing us with a useful reference spectrum. Note also 
that very recent OBE calculations [19| come very close in 
describing the HADES p+p and n+p data consistently. 

The questions now arising are: How does the pair 
yield evolve with increasing system size? Does the in- 
fluence of the hadronic medium set in and, if yes, what 
are its characteristic signals? To address this complex we 
present here results from our measurement of e + e~ pro- 
duction in the medium-heavy reaction system 40 Ar+KCl 
at 1.76^4 GeVin which we also have observed uj produc- 
tion for the first time at SIS18 bombarding energies. In 
section II of this article a brief description of the ex- 
periment, as well as of the analysis procedures, is given. 
In section III the reconstructed e + e~ mass distribution 
is presented and its relevant features in terms of excess 
and vector-meson regions are discussed. In section IV we 
show and discuss other pair observables: transverse-mass 
and angular distributions. Finally, section V summarizes 
and concludes the paper. 



II. THE AR+KCL EXPERIMENT 

The HADES detector operates at the GSI 
Helmholtzzcntrum fur Schwcrionenforschung in Darm- 
stadt with proton and heavy-ion beams being provided 
by the synchrotron SIS18. Technical aspects of the 
spectrometer are described in Q, a schematic view is 
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shown in Fig. [TJ Here we recall that its main component 
serving for electron and positron identification (PID) is 
a ring- imaging Cherenkov detector (RICH). Further PID 
power is provided by (1) the time of flight measured in 
the plastic scintillator time-of- flight wall (TOF), (2) the 
electromagnetic shower characteristics observed in the 
Pre-Shower detector, and (3) the energy loss signals from 
the four wire-chamber planes and the scintillators of the 
TOF wall. A 50 fim thick segmented polycrystalline 
diamond detector (START) placed in front of the target 
provides the precise event time. 




FIG. 1: (Color online) Schematic view of the HADES detector 
as implemented in the simulation. Simulated particle tracks 
are shown as well. 

A beam of 40 Ar ions with a kinetic energy of 
1.756 A GeV was used for the Ar+KCl experiment dis- 
cussed in this paper. The beam intensity was about 
6 x 10 6 particles per 8-second spill. The four-fold seg- 
mented target was made of natural KC1 with a total 
thickness of 5 mm, corresponding to 3.3% nuclear in- 
teraction length. The segmentation of the target (1.3% 
radiation length per segment) helped to minimize the 
conversion of photons into electron pairs in the target 
material. The four segments are indeed nicely visible in 
the reconstructed event vertex distribution along the z 
axis (i.e. beam axis) shown in Fig. [2] 

The on-line event selection was done in two steps: a 
first- level trigger (LVL1) picked out those reactions in 
which the number of hits exceeded 16 in the TOF scintil- 
lators. In total, 2.2 x 10 9 of such events were examined 
with a second-level trigger (LVL2) to find single lepton 
signatures of which about 6 x 10 s were finally recorded. 
In the off-line analysis, events were further filtered by 
cutting on the event vertex reconstructed with a resolu- 
tion of a x — a y ~ 0.4 mm and a z ~ 1.9 mm. These 
vertex cuts removed about 5% of the events in accor- 
dance with the event rate observed in an empty-target 
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FIG. 2: Reconstructed event vertex distribution along the 
beam axis z. The four KC1 target segments are clearly sepa- 
rated. 

run. The LVL1 trigger enhanced the mean pion mul- 
tiplicity approximately two times with respect to the 
minimum-bias multiplicity. From our charged-pion anal- 
ysis and a simulation of the trigger response to UrQMD 
events, we found that this corresponds to a mean num- 
ber of participating nucleons of {A^Yrt 1 ) =38.5 and an 
average charged-pion multiplicity of \{N V + + N^-) = 
3.50 ± 0.12(stat) ± 0.22(sys) (for details see ^ Hj). 

To investigate systematic effects in the dielectron re- 
construction of the HADES experiment, three parallel 
data analyses were done with identification algorithms 
based respectively on (i) a multi-variate analysis (MVA) 
[22l [23| , (ii) a Bayesian approach 0, H3| , and (iii) a 
combination of multi-dimensional selection cuts [3, [25| ■ 
All three analyses led to consistent results with simi- 
lar signal purities. The remaining differences were com- 
pounded with the other systematic uncertainties (see be- 
low) . Identified electrons and positrons were further com- 
bined into pairs. The total number of reconstructed e + e~ 

pairs, N-\ , can be decomposed as iV H = S + CB, 

where S denotes the number of signal pairs and CB 
stands for the number of combinatorial background pairs. 
The former are the correlated e + e~ pairs originating 
from the same parent particle and the latter ones are 
due to combining uncorrelated leptons stemming from 
separate sources, mostly ir° Dalitz decays and external 
photon conversion. To reduce the CB, and hence en- 
hance the S/CB ratio, it is in particular necessary to sup- 
press contributions from photon conversion, from track- 
ing fakes, and from misidentified hadrons. The main 
source of photons are the neutral-pion decays, 7r° — > 77. 
As the induced conversion pairs have mostly small open- 
ing angles, they have been effectively decimated with an 
opening-angle cut requiring a e + e - > 9° in the labora- 
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tory frame. Tracking fakes were suppressed by selecting 
only ring-track combinations of sufficient reconstruction 
and matching quality Q- Additionally, a single- lepton 
momentum cut of 0.1 GeV/c < p e < 1.1 GeV/c confined 
the fiducial acceptance to the region where the combined 
track reconstruction and lepton identification efficiency 
was at least 10%, but typically 30 - 70 %, while the con- 
tamination of the lepton sample by charged pions and 
protons stayed well below 20%. 

Finally, to obtain the e + e~ invariant mass signal, 
the remaining CB was subtracted from all reconstructed 
pairs in the following way: in the low-mass region M ee < 
0.4 GeV/c 2 , where the correlated background from the 
7T° two-photon decay followed by double conversion con- 
tributes most, the combinatorial background was de- 
termined using a method based on like-sign e + e + and 
e~e~ pairs emerging from the same event, i.e. CB — 
2y/N e + e +N e - e -. For lar ger masses, however, where the 
statistics of like-sign pairs is poor, we used a mixed-event 
CB normalized to the like-sign CB [7J. The mixing was 
done between events belonging to the same event class in 
terms of the track multiplicity (five selections) and the 
target segment (four selections), i.e. for a total of twenty 
event classes. This procedure was applied likewise to all 
other pair observables, in particular the pair transverse 
momentum distribution. 
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The resulting invariant mass spectrum of the dielec- 
tron signal, corrected for the detector and reconstruction 
inefficiencies 2 , but not acceptance, is shown in Fig. [3] 
The spectrum is normalized to the avera ge number of 
charged pions - also measured in HADES [2(| - namely 
(N n - + N„+) /2 = 3.5 per LVL1 event. As expected from 
isospin symmetry, this average is a good estimate of the 
actual 7T° yield N v a, i.e. we set N^o = (N v - + N n +) jl. 
The normalization to N^o compensates in first order the 
bias caused by the implicit centrality selection of our 
LVL1 trigger. The spectrum shown represents an av- 
eraged result from the three parallel PID analyses men- 
tioned above. Besides the statistical error bars system- 
atic errors are represented as horizontal ticks. They cover 
systematic effects attributed to the efficiency correction 
and combinatorial background subtraction (20%), to the 
error on the normalization (11%), and to differences be- 
tween the three PID methods (10%). Statistical errors 
are of course point-to-point, the normalization error is 
global, and the remaining systematic errors are slowly 
varying with pair mass. The systematic errors given are 
upper bounds and add quadratically to 25%. 



III. RESULTS FROM AR+KCL 

Here we discuss in more detail the efficiency-corrected 
and CB-subtracted dielectron invariant mass spectrum 
from Ar+KCI (see Fig. [3J. The total yield of ~85k sig- 
nal pairs is distributed over three easily distinguishable 
regions: (i) masses below 0.15 GeV/c 2 , dominated by 
the 7T° Dalitz peak, contribute around 70k, (ii) the in- 
termediate range of 0.15-0.5 GeV/c 2 where the pair ex- 
cess is located, holds about 15k, and (iii) masses above 
0.5 GeV/c 2 where the dileptons from vector meson de- 
cays are expected, total a few hundred pairs only (~ 450). 
All pair observables presented below have been obtained 
from inclusive LVL2-triggered events, i.e. with no further 
centrality cuts. An investigation of different event classes 
selected by analysis cuts on the hit multiplicity revealed 
indeed a slight dependence of the normalized pair yields 
on centrality. However, as in this still rather small reac- 
tion system only a limited range of A part can be scanned 
via such multiplicity cuts, we discuss below the A part de- 
pendence only in the context of a comparison of Ar+KCI 
with C+C. 



FIG. 3: (Color online) Reconstructed e + e~ mass distribu- 
tion in Ar+KCI collisions (averaged over three PID analy- 
ses, efficiency-corrected, CB subtracted, and normalized to 
N^o). Statistical and systematic errors of the measurement 
are shown as vertical bars and horizontal cups, respectively. 
Curves represent the 7r° and rj Dalitz components, as well 
as the u) contribution (Dalitz and direct) simulated with the 
event generator Pluto. Also shown are the excess yield over 
the simulated cocktail (shaded area) and a fit (exponential 
+ Gaussian curves) to the data in the mass range 0.25 - 
0.9 GeV/c 2 (see III.C for details). 



A. Low-mass pairs 

The low-mass region contains the bulk of the pair 
yield, but it is also the one most strongly affected by the 



2 Inefficiencies were determined with an event overlay technique in 
which simulated lepton tracks were embedded into real events, 
reconstructed, and tallied. 
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momentum-dependent efficiency corrections. As more 
than 90% of this yield stems from the 7r°, it offers a conve- 
nient handle for validating our dielectron reconstruction 
and normalization procedures. To do this check we simu- 
lated the pion and eta contributions to the e + e~ cocktail 
with the Pluto event generator [2(| H3] using the meson 
multiplicities and source parameters given in Tab. HI In 
case of the ir° they were taken from our own charged- 
pion data 0, HJ, in case of the 77 they were interpo- 
lated from TAPS two-photon measurements of 1.5 and 
2.0 A GcV Ar+Ca (Ca+Ca) reactions [HG3. Note that 
for both mesons a two-slope parameterization has been 
used. The validity of this interpolation is confirmed in 
Fig. U where ir° and rj mid-rapidity l/m\dr N/dm ±dy 
spectra from TAPS are compared with the correspond- 
ing 7r + and 7r~ average measured in Ar+KCl (2fjj (mid- 
rapidity y — 0.858 and \y — yo\ < 0.1). In this figure, 
the different centrality selections of the TAPS (minimum 
bias) and HADES (LVL1) experiments are compensated 
by normalizing to the mean number of participating nu- 
cleons (A par t) = 20 and 38.5, respectively. The Ar+KCl 
charged-pion average falls nicely in between the neutral 
pion (and eta) data, as expected from the smooth beam 
energy dependence of pion production. The characteris- 
tics of the vector meson sources (p, w, <j)) are still largely 
unknown and have tentatively been set as given in Tab. [I] 



TABLE I: Thermal source parameters used in the Pluto sim- 
ulation of the dielectron cocktail. For all listed particles we 
give the multiplicity (N) within our LVL1 trigger condition, 
the source temperatures (Ti and Tb), the relative strength 
(frac = ci/ (ci +C2), where ci and C2 are the amplitudes of the 
two components) of the first component, the polar anisotropy 
(A2), and the helicity coefficient (B) of the dielectron decay 
(see Sec. III.E). 



Part. 



7T 
V 

A+-° 

ui 

P 



TV Ti [MeV] T 2 [MeV] frac A 2 B 



3.5 
8.8 ■ 1(T 2 

3 TV(tt°) 
6.7- 1(T 3 

5 • 1(T 2 
2.6 • 1(T 4 



52 
52 
80 
80 
80 
80 



89 



0.85 
0.85 

1 

1 

1 

1 



0.75 1 

0.75 1 

0.75 1 

0.75 

0.75 

0.75 



Turning back to Fig. [3j one can see that the low-mass 
pair yield is indeed described very well by our calculation, 
lending confidence to the reconstruction process. 



B. Intermediate- mass excess 

Comparing a Pluto simulation of long-lived sources 
(i.e. emitting mostly after freeze-out) with the data in 
Fig.[3]reveals that also in Ar+KCl the contributions from 
r] (and uj) Dalitz decays do not exhaust the measured pair 
yield at intermediate masses, i.e. for M ee ~ 0.15 — 0.5 
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FIG. 4: (Color online) Average of the mid-rapidity charged 
pion d 2 N/dm±dy distributions, N = l/2[A r 7r+ + N„~], mea- 
sured with HADES in the 1.76 A GeV Ar+KCl reaction 
[20I [2lT ] (squares linked by solid curve), compared with the 
corresponding neutral pion and eta data from 1.5A GeV 
40 Ar+ nat Ca and 2.0+ GeV 40 Ca+ nat Ca reactions obtained 
with the photon calorimeter TAPS [ll|, 03] (circles and tri- 
angles). All yields are normalized to their respective {A par t)', 
error bars shown are statistical. 



GeV/c 2 . Just like in our previous C+C measurements 
[1, Isj , there is a strong excess over the cocktail of known 
long-lived sources. We know furthermore from our com- 
parison [l8| of dielectron production in free nucleon- 
nucleon and C+C reactions that in the light carbon- 
carbon system the excess yield N exc is of baryonic origin: 
A decays and TV TV - i.e. mostly pn - bremsstrahlung. The 
C+C reaction seems to be in first order an incoherent 
superposition of elementary NN processes. In addition, 
whereas between 1 and 2A GeV rj production increases 
steeply with bombarding energy (from sub-threshold to 
above threshold), the excess yield rises like pion pro- 
duction, i.e. only mildly. This means that our isospin- 
averaged ^[pp + np] pair spectrum measured at 1.25 GeV 
[HI can serve as TV TV reference for the ^-subtracted pair 
yield observed in the present 1.76 A GeV Ar+KCl run as 
well, both normalized to their respective ir° multiplicity. 
Note, however, that this reference is of course available 
only up to its kinematic cutoff at M ee = 0.55 GeV/c 2 , 
corresponding to the 1.25 GeV bombarding energy used 
in the TV TV experiments. 

Figure[5](top) shows the Ar+KCl e + e~ invariant-mass 
distribution after subtracting the simulated 77 component 
and normalizing to N^o , together with the TV TV reference 
from [lH (adjusted to the acceptance, i.e. magnetic field 
and momentum cuts, of the present experiment), also rj 
subtracted and normalized to its ir° multiplicity (aver- 
aged from p+p and p+n data). In this comparison we 
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FIG. 5: (Color online) (a) Comparison of the Ar+KCI 
invariant-mass distribution with an isospin-averaged reference 
from p+p and n+p data [l8| . For clarity systematic error 
bars are shown only on every second data point (vertical bars 
are statistical, cups are systematic). Both data sets are nor- 
malized to their respective pion multiplicity and have their 
respective 77 Dalitz yield subtracted. The dashed lines are 
meant to guide the eye. (b) Ratio of the heavy-ion mass dis- 
tributions (Ar+KCI and C+C) to the 1/2 [pp+np] reference, 
whose total error (statistical and systematic added quadrati- 
cally) is indicated by the shaded band. Note that all data sets 
are shown within the acceptance of the Ar+KCI experiment. 



the yield from the heavy system exceeds the NN refer- 
ence by a factor of ~ 2.5 — 3. This is also visible in the 
lower part of Fig. [5] where the ratios of the following pair 
yields are shown: Ar+KCI/ N+N, and C+C /N+N for 
1 and 2A GeV. For this the C+C data were taken from 
[H HI and transformed into the acceptance of the present 
experiment. The Ar+KCI/ N+N ratio is very close to 
unity at low masses, dominated by the 7T° Dalitz pairs, 
but for M > 0.15 GeV/c 2 it rises to about 3, indicating 
the onset of processes not accounted for in the reference 
system. Both representations prove that a qualitative 
change happens in the nature of the excess yield when 
going to the heavier system. Consequently, in contrast 
to the C+C system, Ar+KCI can not anymore be seen 
as a superposition of independent N+N collisions. A 
more complex picture involving multi-body and multi- 
step processes and maybe even in-medium modifications 
of the involved hadrons is required. Note also that a 
scaling with the number of binary nucleon-nucleon col- 
lisions N co u might be more appropriate to describe the 
observed variation of the excess yield with system size. 
Indeed, (N co u) calculated within a Glauber approach [28[ 
increases faster than (A part ) when going from our LVL1 
C+C to LVL1 Ar+KCI events, namely by a factor 6.1 for 
(N coU ) vs. 4.5 for (A part ). 

Combining the dielectron results from HADES and 
from the former DLS experiment we can now study 
the evolution of the excess over cocktail with beam en- 
ergy and system size. To do so we have compiled in 
Fig. |5] the excess yields integrated over the mass region 
M ee = 0.15 — 0.5 GeV/c 2 from all available reaction 
systems 0, H 0- For comparison, inclusive n° and r\ 
multiplicities measured in photon calorimetry with the 
TAPS detector [H El are plotted as well. Note that all 
yields are extrapolated to the full solid angle 3 and are 
normalized to their respective average A part in order to 
compensate for differences in the centrality selection of 
the various experiments. The normalization also takes 
out the trivial system-size dependence of the yields, as 
visible from the closeness of the C+C and Ca+Ca me- 
son curves . The somewhat smaller pion multiplicity per 
Apart of Ca+Ca can be attributed to meson re-absorption 
in this larger system. Note, however, that the eta multi- 
plicities start out with the opposite behavior at low beam 
energy and switch only around Eb earn =1.5 A GeV to the 
absorption-dominated scaling. This crossing can be ex- 
plained by the transition from the sub-threshold regime, 
where multi-step processes favored by a larger reaction 
volume are important [29[, to above threshold produc- 
tion. 

Next one can see that the dielectron excess follows pion 
production with rising bombarding energy, as we stated 



do not correct for the slight isospin asymmetry of the 
Ar+KCI system (N/Z = 1.15). Due to the normaliza- 
tion and the uae of a common acceptance both distribu- 
tions agree in the 7r° Dalitz peak. They differ, however, 
strikingly for masses between 0.15 and 0.5 GeV/c 2 where 



3 Assuming similar geometric acceptances for excess pairs and r\ 
Dalitz pairs. 

4 We consider here the systems Ar+KCI and Ca+Ca as being 
equivalent in size and isospin. 



7 



CO 
CL 



10" : 



10" 



10" 



7T 




Ca+Ca 



6 _ 



".Zr" 
e*e excess 

0.15 <M ee < 0.5 ^ c+C 



Ar+KCI 



*-■ 







1 



"beam 



[A GeV] 



FIG. 6: (Color online) Inclusive multiplicity per participant, 
N 4n /(Ap ar t), as function of beam energy of the e + e~~ pair 
excess over the r\ Dalitz yield, and of ty° and r\ production in 
heavy-ion reactions. The excess yield, defined in the mass 
range M ee = 0.15 — 0.5 GeV/c 2 and extrapolated to 4-7T, is 
shown for HADES Ar+KCI and C+C data (full triangles) 
[3, @] as well as for DLS data (open triangles) [f| . The n° 
and r) results are from TAPS measurements in C+C (squares, 
solid curves) and Ca+Ca (circles, long-dashed curves) colli- 
sions [lot Q2] ■ The curves are polynomial fits to these data 
used to interpolate the multiplicities as a function of bom- 
barding energy (see [11J). For easier visual comparison with 
the energy dependence of the dielectron excess the latter is 
overlayed with the 7r° curves (short-dashed) down-scaled by 
factors of 1.8 ■ 10" 5 for C+C and 4.3 ■ 10" 5 for Ar+KCI. 



already before [8|. This turns out to be true for both 
the C+C and Ca+Ca collision systems, as one can see 
from the excellent match with the scaled-down pion pro- 
duction curves in the figure. Such a behavior has been 
interpreted as being characteristic of a production mech- 
anism not driven by the excitation of heavy resonances, 
but rather by low-energy processes like pion production 
and propagation involving the A and, maybe, low-mass 
p excitations as well as bremsstrahlung jl8| . 

As already pointed out, the systematics of excess yields 
has become sufficiently rich to allow also for a study of the 
system-size dependence of the electromagnetic radiation 
from the nuclear medium. Above we concluded that the 
comparison of the excess yield obtained in Ar+KCI with 
the NN reference reveals a non-trivial behavior of pair 
production. This is also supported by Fig. [6] where an 
increase of a factor ~ 2 is visible when moving from the 
C+C to the Ca+Ca system. Evidently the excess yield 
must scale faster than linear with A part in contrast to the 
behavior of e.g. pion production in heavy-ion reactions. 
The mass dependence can indeed be further quantified 
by adjusting a N exc oc A part scaling law to the yields. 
Using the Ar+KCI excess obtained at 1.7&A GeV and in- 



terpolating the C+C excess for that beam energy we get 
a coefficient a — 1.41 jig 27- Using instead the DLS point 
measured at 1.04^4 GeV results in a similar scaling co- 
efficient. Note that, when varying (A part ) by comparing 
systems of different size, the corresponding scaling expo- 
nent for pion production has been found to be a ~ 0.85 
[l2l [30l |. independent of beam energy, and the one for 
eta production to be a ~ 1.2 at 1A GeV, decreasing to 
a ~ 0.8 at 2^4 GeV [l2j]. These numbers confirm that the 
dielectron excess scales with system size very differently 
than the freeze-out yields of pions and eta mesons. 

All of our observations put together may be interpreted 
as the onset multi-body and multi-step processes in the 
hot and dense phase created in collisions of nuclei of such 
a size. The penetrating nature of the dilepton probe of- 
fers then a natural explanation for the behavior of the ex- 
cess e + e~ yield if one keeps in mind that for a sufficiently 
large number of participating nucleons or, in other terms, 
for a sufficiently large collision volume the detected ra- 
diation is the integral over the full time of the complex 
heavy-ion reaction and not just a snapshot at freeze-out. 
It will be interesting to follow up on this trend when still 
heavier systems are added to the systematics. 



C. High- mass pairs from omega decays 
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FIG. 7: (Color online) Zoom into the measured e + e~ mass 
distribution in Ar+KCI (efficiency corrected and normalized), 
together with a least-squares fit of the 10 vector meson peak 
by a sum of a Gauss function and an exponential background 
(see text for details). Error bars are like in Fig. [3] 

We turn now to the high-mass region of the invariant 
mass spectrum and, in particular, to the clearly protrud- 
ing peak structure which we attribute to the direct decay, 



8 



oj — > e + e~ , of the omega vector meson. A linear zoom- 
in onto the vector-meson region is shown in Fig. [7] The 
peak visible at the oj pole position holds about 40 recon- 
structed pairs, limiting unfortunately the extent to which 
one can possibly go in its analysis. These data constitute 
the very first observation of omega production in a heavy- 
ion reaction at such a low beam energy, in fact, an energy 
even below the production threshold in free N+N colli- 
sions (E$£ = 1.89 A GeV). One expects that most of 
the omegas contributing to this peak decayed after hav- 
ing left the reaction zone, i.e. after freeze-out. Recently 
measured u> photoproduction cross sections [3lT[33| have 
been interpreted [l| in the sense of a strong broadening 
(up to 150 MeV) of the decay width of this meson in 
the nuclear medium already at normal nuclear density. 
We do not observe such a modification in our omega sig- 
nal: the shape of the observed peak is solely determined 
by the detector response, i.e. by the intrinsic momen- 
tum resolution of the HADES tracking system. In this 
mass region also p° decays and baryonic resonance de- 
cays are expected to contribute to the dielectron yield, 
but they add up to a broad continuum underneath the 
omega peak. For masses above 0.9 GeV/c 2 the statistics 
is running out quickly and there is no recognizable struc- 
ture at the pole position of the 4> meson (M^ = 1.019 
GeV/c 2 ). 

All of this justifies fitting the whole mass region with 
the sum of a Gauss shape and an exponential function, as 
shown in Fig. [7] The fit (x 2 /ndf = 11.8/18) provides a 
peak position of M u = 0.770 ± 0.011 GeV/c 2 , a width of 
a u = 0.022±0.010 GeV/c 2 , and an integrated signal over 
the continuum corresponding to (3.9 ± 1.7) • 10 -8 . The 
peak centroid agrees hence within about one standard de- 
viation with the listed oj pole position at 0.783 GeV/c 2 
[34j . Furthermore, detector simulations show that part 
of the observed down-shift (~ 10 MeV) is due to the 
combined energy loss of the electron and positron in the 
inner part of the HADES detector. The peak width is 
dominated by the HADES mass resolution a/M at the oj 
pole mass of 3%. Finally, its integral has been corrected 
for the branching ratio of the direct e + e" decay [34j as 
well as for the acceptance of 0.29 (obtained from a Pluto 
simulation done for a thermal source with a temperature 
of T = 84 ± 2 MeV, as found in our <p -» K+K~ analy- 
sis [35[ for the 4> meson). This resulted in a normalized 
yield of N^/N^o = (1.9 ±0.8) • 10~ 3 , corresponding to an 
omega LVL1 multiplicity of M^ VL1 = (6.7 ± 2.8) • 10~ 3 . 
Fits with more sophisticated peak shapes taking into ac- 
count the slightly asymmetric momentum response of the 
detector gave very similar results. The acceptance cor- 
rection depends mildly on the phase-space distribution 
used in the Pluto simulation: it ranges from 0.34 at 
T = 50 MeV to 0.24 at T = 140 MeV (see also the 
discussion of the pair m± slopes in the next subsection) . 
It depends even less on the assumed polar distribution: 
5% decrease when varying Ai from to 1. All those ef- 
fects are finally subsumed into an additional systematic 
error on the multiplicity of 25%. With the oj yield known, 



both its contributions - Dalitz and direct - to the pair 
cocktail can be simulated in Pluto; they are shown to- 
gether with the mass spectrum in Fig. [3] The oj decays 
contribute evidently only a small part to the total pair 
yield at intermediate and low masses. Note finally that 
the average oj momentum in the nucleus-nucleus center- 
of-mass within the HADES acceptance is found from our 
data to be p = 0.43 GeV/c. This is at least a factor 
two smaller than the momenta t ypi cally observed in oj 
photoproduction experiments [3ll433j . 
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FIG. 8: (Color online) Comparison of the R^/u, ratio obtained 
in this work with its statistical model (THERMUS fit) value 
as well as with a compilation of data from elementary p+p and 
7T+N reactions (see text). The ratio is plotted as a function 
of the excess energy e in the NN — > NN<^ and the 7rN— >N0 
reactions, respectively. 

The oj multiplicity can be discussed in the context of 
either a scenario of complete thermalization at freeze- 
out or, in the other extreme, of production in elemen- 
tary N+N collisions. As HADES is a general-purpose 
charged-particle detector, besides the dielectron results 
presented here, a wealth of information has been obtained 
as well on hadron production in Ar+KCI. These findings 
have already been published in I20l on n^, in [35[ on 
K + , K~ , and </>, in [21( on K®, in [36| on S - , and finally 
in [33 on A and S ± . 

In particular, from our K + — K~ correlation analysis 
35], a LVL1 multiplicity of M^ = (2.6 ± 0.7(stat) ± 
0.1 (sys)) ■ 10~ 4 has been found as well as a transverse- 
mass slope at mid-rapidity of = 84 ± 8 MeV. To- 
gether with the oj multiplicity, this gives a <p/oj ratio 
of = 0M3t O o o Q l° 5 (stat) ± O.Oll(sys). The exper- 

imental ratio can be compared to various predictions, 
running from pure m± scaling in 47r solid angle, giv- 
ing R ~ 0.042, to a full-fledged statistical hadroniza- 
tion model calculation performed with the THERMUS 
code [l^ fitted to our hadron yields [37} and resulting in 
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R = 0.063 ± 0.008. Hence, statistical descriptions agree 
within error bars with the experimental i u . As already 
discussed in ref. j^, the THERMUS model does well in 
reproducing our measured hadron yields, including those 
of particles with open or hidden strangeness, with the 
notable exception of the double-strange S~ which, how- 
ever, at 1.76A GeV is produced far below its threshold 
of 3.57 GeV in free NN collisions. 

The opposite extreme to complete thermalization is 
given by elementary N+N and n + N reactions where 
the 4>/uj ratio is traditionally investigated in the context 
of the so-called OZI rule violation [39-41]. The ratio 
obtained in those reactions for small values of the excess 
energy (e = E c . m . — Ethr) notoriously exceeds predictions 
based on the <f> — lu mixing angle and is sometimes related 
to a possible ss admixture in the nucleon ground-state 
wave function. Figure |8] shows R^/u obtained in this 
work and the THERMUS value from a fit to HADES 
data together with results from elementary p+p [HI, |43[ 
and 7r + N [44[ reactions, all plotted as function of the 
excess energy in the NN— »NN(/> and 7rN— >Ncf> reactions, 
respectively. This is different from the common definition 
in literature where the <fr and w yields are both taken at 
the same excess energy, corresponding hence to different 
bombarding energies, whereas we take the ratio of yields 
measured at a common beam energy. In fact, to do this, 
we divided the measured <f> cross sections by an interpo- 
lation of the omega cross sections based on the parame- 
terization proposed in [40| and, in case of the p+p data, 
updated in [45j . One can see from the comparison that in 
the heavy- ion reaction the ratio R^ / u is more than an or- 
der of magnitude larger than in NN collisions and also at 
least a factor 3-5 larger than in pion-induced processes. 
One should furthermore keep in mind that mostly low- 
momentum pions are produced in 1 - 2A GeV heavy-ion 
reactions while the <f> production threshold is at — 
1.56 GeV/c; in N+N collisions the production threshold 
is at 2.60 GeV. Consequently, the <f> meson is produced 
sub-threshold here (e < 0) and more complex, multi-step 
processes involving short-lived resonances 46] and/or hy- 
perons [47| might contribute. On the other hand, the 
ratio could also be influenced by final-state effects of the 
vector mesons in the nuclear medium (see e.g. [l] for a 
discussion of u> and <f> absorption). In the end, our obser- 
vation seems to support the picture of meson production 
in a rather long-lived and thcrmalizcd fireball. 



In all other regions, however, multiple sources contribute 
and their disentanglement is not trivial. We have for- 
tunately constraints on the r\ Dalitz yield from earlier 
TAPS measurements and now also on the uj Dalitz yield 
from our own analysis of the multiplicity of this particle 
(see discussion in III.C). 
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D. Dielectron m± distributions 

We present now phase-space distributions of e + e~ 
pairs in Ar+KCl. When discussing the pair mass spec- 
trum (see Fig. [3]) we distinguished different mass regions 
of interest. Indeed, the pair spectrum is a complicated 
cocktail emitted from various processes and at different 
phases of the heavy-ion reaction. As pointed out above, 
at low masses the situation is rather simple because this 
region is dominated by pairs from the 7r° Dalitz decay. 



FIG. 9: (Color online) Reconstructed rapidity (a) and trans- 
verse mass at mid-rapidity (b) distributions of e + e~ pairs 
with M ee < 0.13 GeV/c 2 . Data are normalized to iV^o as 
well as corrected with the detector efficiency and acceptance. 
The error bars are statistical. The histograms correspond to 
a simulated (Pluto) cocktail of thermal sources. A Gauss 
fit to dN/dy (top, solid curve), shown as well, results in 
(y) = 0.83 ± 0.03 and a y = 0.91 ± 0.07. 

To characterize the dielectron yield beyond its mass 
distribution one has to reconstruct other pair observ- 
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ables, in particular its longitudinal and transverse phase- 
space population. The longitudinal dimension is usually 
covered by plotting rapidity density as function of ra- 
pidity, dN/dy, and the transverse one by plotting either 
dN/dp± or dN/dm± with m± = \fp\ + M| e as function 
of p± or m± . The thermal nature of a particle source can 
be best recognized by plotting either its 1 / m\ dN/ dm± 

distribution at y = yo or its 1 / m ! 2 dN/ dm± distribu- 
tion integrated over all rapidities [481 ] . Indeed, in semi- 
logarithmic representation and for m± >> T, both func- 
tions turn into a straight line where the inverse-slope pa- 
rameter T may be interpreted as the source temperature. 
In case of measuring dilepton pairs the situation is further 
complicated by the decay kinematics of the three-body 
Dalitz decays. Consequently, the exact nature of the par- 
ent distribution can be distorted in the observed e + e~ 
distribution. Note that, in order to obtain meaningful 
slopes, these distributions have to be corrected not only 
for efficiency but also for acceptance including the detec- 
tor geometry as well as momentum and opening angle 
cuts. As mentioned in the discussion of the uj multiplic- 
ity determination, the acceptance correction has been ob- 
tained from Pluto simulations of a full pair cocktail (with 
the source parameters listed in Tab. H]), while varying its 
source parameters to quantify systematic effects. The 
pair acceptance has thereby been determined as a one- 
dimensional function of transverse mass, averaged over 
rapidity within a given mass bin, and vice- versa. We have 
verified that this procedure gives results compatible with 
the more complex multi-dimensional correction as func- 
tion of mass, transverse momentum, and rapidity. The 
resulting normalized dN/dy and dN/dm± spectra are 
shown in Fig.[S]for pairs of M ee < 0.13 GeV/c 2 , together 
with the corresponding simulated spectra. This low-mass 
bin - being dominated by ir° Dalitz yield - is described 
to better than 10% by our Pluto event generator, as seen 
from the overlayed histograms. This agreement further 
strengthens our confidence in the pion phase-space distri- 
bution used in the simulation as well as in our dielectron 
reconstruction procedures in general. 

In the context of this analysis we have also done a care- 
ful investigation of the signal purity, as one might fear 
that particularly the high m± pairs could be contam- 
inated by misidentified high-momentum hadron tracks 
and/or fake tracks. This purity study has been done with 
an event mixing technique and confirmed as well with full 
simulations of the reconstruction and particle identifica- 
tion. Indeed, while our lepton purity is on average better 
than 0.95, it decreases with increasing lepton momen- 
tum, resulting nonetheless in a dielectron purity which 
remains better than 0.7 up to m± values of 1.5 GeV/c 2 . 
Note that hadron and fake impurities in the lepton sam- 
ple lead to uncorrelated pairs only, and thus increase the 
combinatorial background which is of course subtracted, 
as discussed in section II. We have checked in simulations 
that the CB subtraction indeed removes these additional 
uncorrelated contributions. 

To take advantage of our full pair statistics, we have 
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FIG. 10: (Color online) Reconstructed pair l/rn 3 / 2 dN/drn ± 
distributions, normalized to the 7r° multiplicity, for the full 
rapidity range and different mass selections given in the l.h.s. 
legends (in MeV). Efficiency and acceptance corrections are 
applied; error bars are statistical. Exponential fits to the 
high-mj_ region of the data are shown as dashed curves with 
the corresponding inverse-slope parameter given (in MeV) in 
the second line of the legends. Note also the scaling factors 
(in parentheses). 



opted to use the 1/m 3 / 2 dN/dm± = N exp(—m±/T) 
representation in our systematic investigation of the 
transverse momentum distribution for several bins of 
pair mass displayed in Fig. 1101 With increasing pair 
masses contributions from r\ Dalitz, A (and N*) Dalitz, 
bremsstrahlung, and finally u), p° and (very few) <f> decays 
are successively probed. Although the limited statistics 
of our data required rather wide bins, particularly for the 
highest masses, one can see a distinctive pattern emerge: 
as one progresses from low to higher M ee , the slope of the 
pair transverse-mass spectra first remains approximately 
constant at T around 80 MeV, but when approaching the 
vector meson region, it rises steeply to reach a value as 
high as 130 MeV. 

While the first mass bin is dominated by ir° Dalitz 
pairs, as emphasized in discussing Fig. [9j the next two 
bins cover the intermediate- mass region (0.15 < M ee < 
0.5 GeV/c 2 ), with contributions from 77 Dalitz, A Dalitz, 
NN bremsstrahlung and maybe other sources. This is 
the region of the pair excess which we would like to char- 
acterize as much as possible. To do this, we have again 
subtracted the eta component simulated with Pluto by 
making use of the known r\ multiplicity and source tem- 
perature. The resulting excess dN/dm± distribution is 
shown in Fig. [TT] together with corresponding data ob- 
tained in C+C at 1 and 2A GeV, as well as with a 
reference from elementary nucleon-nucleon collisions, ob- 
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FIG. 11: (Color online) Reconstructed l/m ± dN/dm± distributions of pairs with 0.15 < M ee < 0.5 GeV/c and — oo < y < oo 
in C+C at 1A GeV (a), in Ar+KCI at 1.76A GeV (b), and in C+C at 2A GeV (c). Reference spectra from elementary NN 
collisions are also shown, namely the average of p+p and n+p at 1.25 GeV (open crosses), and p+p at 2.2 GeV (open triangles). 
Efficiency and acceptance corrections are applied. All distributions have their respective r) contribution subtracted and are 
normalized to their respective pion multiplicity N n o . Error bars are statistical. Dashed lines are exponential fits, with the 
corresponding inverse slope parameter given in MeV. 



taincd from the average of our p+p and n+p results at 
1.25 GeV [l|[ as discussed in III.B. The spectrum from 
elementary p+p collisions at 2.2 GeV [49j is also shown, 
but at this energy, unfortunately, the corresponding n+p 
yields needed for isospin averaging are not available. All 
distributions are normalized to their respective neutral 
pion multiplicity, N w o , and have their respective i] Dalitz 
contribution subtracted 5 . The figure shows that - within 
error bars and up to m± ~ 0.5, respectively ~ 0.8 - the 
light C+C system behaves at both bombarding energies 
very much like the NN reference: normalized yields and 
slopes agree to a large extent. Note again that the ref- 
erence spectrum can have yield only up to its kinematic 
cutoff at m_L = 0.55 GeV/c 2 (0.89 GeV/c 2 for 2.2 GeV 
p+p). In contrast to C+C, however, the Ar+KCI system 
displays a large excess over the elementary reference, just 
as found already in the comparison of the pair invariant 
mass spectra. 

Moving finally to the large kinetic slopes found in the 
two upper mass bins, we note that this observation is sur- 
prising and difficult to reconcile with the assumption of 
a completely thermalized particle source. As one moves 
away from the three-body decays dominant below the 77 
Dalitz edge at 0.547 GeV/c 2 , the two-body decays of the 
vector mesons contribute more and more, and one ex- 
pects indeed an increase of the slope parameter. This is 
just a natural consequence of the decay kinematics. On 
the other hand, from our THERMUS fit to the full set 
of hadron yields measured in Ar+KCI, we have found a 
chemical freeze-out temperature of the fireball of 76 MeV 



mpe : 

(see section III.C and [33]) and kinetic slope parameters 



5 For 2.2 GeV p+p only the exclusive 77 production has been sub- 
tracted. 



in the range of 70 - 95 MeV [33], i.e. of similar magnitude. 
In particular, the slope at mid-rapidity of the <fi meson 
was found to be 84 ± 8 MeV in the K + K~ channel (35j . 
Unfortunately our statistics is not large enough to allow 
a tight selection around the w pole mass and thus obtain 
the slope for a clean uj sample. From transport calcula- 
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FIG. 12: (Color online) Comparison of dielectron 
1/m 3 / 2 dN/dm± distributions with M ee > 0.5 GeV/c 2 and 
—00 < y < 00 in Ar+KCI and C+C (see text for details). 
Efficiency and acceptance corrections are applied. The distri- 
butions are normalized to the respective reaction pion yield 
N^o. Error bars are statistical. Dashed lines are exponential 
fits, with their corresponding inverse slope parameter given 
in MeV. Note also the scaling factors. 
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tions [50|| we estimate that the pair cocktail selected by 
our uppermost mass window (0.65 - 1.2 GeV/c 2 ) contains 
sizeable contributions from the vector mesons p, uj and 4>, 
but its true composition remains of course uncertain. 

Analyzing furthermore the corresponding m± slopes 
from the data sets of the C+C system Q , and compar- 
ing them with the present Ar+KCl result (see Fig. IT2|) , 
we find a comparatively large slope in the light C+C sys- 
tem at 2A GeV. At the lower beam energy of 1A GeV, 
however, the C+C slope is found to be much smaller, 
but still large when compared with the 58 MeV slope of 
charged pions observed in this system [5l|. 

Presently we can only speculate about various effects 
that can lead to such a behavior of the transverse mass 
distributions. For example, collective effects, like radial 
flow, produce large effective temperatures and, in fact, 
yield slopes increasing with particle mass. However, such 
a trend is not visible in the systematics of kinetic slopes 
that we observed in Ar+KCl [37j and, moreover, in this 
rather small s yste m the radial flow is not expected to 
be important [52j . Another mechanism that has been 
proposed is linked to the final-state interactions of the 
produced vector mesons. Indeed, if these interactions 
are strongly momentum-dependent, they can modify the 
spectral distributions. This has also been proposed as an 
explanation for the depletion of uj yield observed at low 
p_i_ in In+In collisions by the NA60 experiment at the 
CERN-SPS [53j]. Re-absorption cross sections of the uj 
meson in cold nuclear matter have been calculated in a 
OBE approach by the authors of ref. |54[ and they found 
them indeed large and strongly momentum-dependent. 
Inserting these cross sections into transport calculations, 
they also predicted observable effects in the transverse- 
mass spectra of uj produced in heavy-ion collisions. Note 
that re-absorption of the vector mesons is closely related 
to their collisional broadening in the nuclear medium. 
From recent measurements of the transparency ratios 
in uj [H, [55[ and (j> [HH] photoproduction the collisional 
broadening of both mesons has been found to be quite 
large (at p — po, = 100—150 MeV), although its exact 
momentum-dependence could not yet be sufficiently con- 
strained (see [l| for a discussion). Yet another factor that 
might influence the shape of the pair transverse mass dis- 
tribution are the spectral functions of the various dilep- 
ton sources contributing. In particular, any enhancement 
at large masses due to increased form factors, as pre- 
dicted e.g. for resonance decays within vector-dominance 
models [57], HH , could lead to spectral distortions tran- 
scended in their characteristic m± slopes. To conclude, 
the interpretation of the large dielectron slope parame- 
ters observed in our heavy-ion data remains challenging. 



E. Dielectron angular distributions 

Angular distributions of the emitted dielectrons consti- 
tute yet another observable of interest. Various emission 
angles can be reconstructed. Here we focus on two spe- 



cific ones: (i) the center-of-mass polar angle c . m ., i.e. the 
angle between the direction of the virtual photon in the 
A + A reference frame and the beam axis, and (ii) the so- 
called helicity angle a, i.e. the angle between the direction 
of the virtual photon in the reference frame of the mother 
particle (e.g. 7r°,ry, A for the three-body Dalitz decays 
and the fireball for the two-body direct decays) and the 
direction of the electron (or positron) in the pair frame. 
This particular choice of the latter angle corresponds to 
its definition in the so-called Jacob- Wick frame [HI, [6(| ■ 
Technically it requires a double Lorentz transformation of 
the lepton momenta: first from the laboratory frame into 
the parent particle frame and second from there into the 
virtual photon frame. Considering the decay kinematics, 
one can convince oneself that the polar angle of the dielec- 
tron will have at least reminiscence of the polar emission 
angle of the mother particle for three-body decays (i.e. 
Dalitz decays) and be equal to it in case of a two-body 
decay (i.e. vector meson — > e + e~). Likewise, the recon- 
struction of the helicity angle is exact only for two-body 
decays and approximate for three-body decays, because 
the third product goes undetected in the inclusive e + e~ 
reconstruction we did. Nonetheless, our simulations show 
that these angular distributions are not completely atten- 
uated if one makes the approximation that the decaying 
particle is at rest in the nucleus- nucleus center-of-mass 
frame and thus information can be gained on the parent 
particle and its decay. The amount of attenuation de- 
pends on the specific decay kinematics and on the source 
temperature of the mother particle. 

The helicity distribution is of particular interest as it 
allows to probe the degree of polarization of the virtual 
photon. It can be proven that pseudoscalar Dalitz decays 
are self-polarizing [gj [UJ and lead to helicity distribu- 
tions of the form dN/da cx 1 + £>cos 2 a with B = 1, 
where a is the helicity angle. This expectation has been 
confirmed long ago for the 7T° in a study using the charge- 
exchange reaction n~p — > Tr°n [63| and again more re- 
cently in exclusive p+p measurements at 2.2 GeV per- 
formed by the HADES collaboration [64]. This HADES 
measurement provided in addition the very first observa- 
tion of the helicity distribution in 77 Dalitz decays. The 
simulations we have done to determine the sensitivity of 
our heavy-ion data to these effects reveal that the at- 
tenuation, caused by the incomplete reconstruction of 
the Dalitz decays, reduces the helicity anisotropy coef- 
ficient from unity to B ~ 0.6 — 0.7 for ir° — > 7e + e~ and 
77 — > 7e + e~ decays, i.e. leaving it still quite sizeable. In 
case of the A —> Ne + e~ decay, the calculation of the he- 
licity distribution is much more involved, but it has been 
done and likewise a distribution of type 1 + B cos 2 a with 
B ~ 1 is expected [65j]. Exclusive dielectron data taken 
with HADES in 1.25 GeV p+p collisions have also con- 
firmed the latter prediction [66| . Finally, in the two-body 
decays - p°,uj,<fr — > e + e~ - the dilepton simply carries 
the full polarization of the vector meson. Measuring he- 
licity angles might hence help to unravel the different 
components of the pair cocktail [65], |67T - [69| . These ideas 
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FIG. 13: (Color online) Acceptance of the HADES detector 
for the dielectron nucleus-nucleus center-of-mass polar emis- 
sion angle 9 c . m . (a) and helicity angle a (b), represented for 
the three mass selections indicated in the figure. 



have also been applied recently by the NA60 experiment 
to characterize the thermal nature of high-mass dimuon 
radiation emitted in high-energy heavy-ion reactions at 
the CERN-SPS [z|. 

Experimentally the dielectron angular distributions are 
only obtained within the detector acceptance shown in 
Fig. [T21 which they need to be corrected for. We have 
done this by dividing a given reconstructed polar distri- 
bution with a corresponding simulated cocktail distribu- 
tion for which isotropic emission of the parent particle 
was assumed. Likewise, the experimental helicity dis- 
tributions were divided by the corresponding simulated 
cocktail distribution assuming an isotropic emission of 
the decay lepton. From our simulations we expect that 
such a ratio, besides correcting for the acceptance, will 
reveal deviations of the data from isotropy. The ir° domi- 
nated low-mass pairs can again serve as a test bed for the 
procedure. For these one expects to observe a polar dis- 
tribution reminiscent of the known pion polar anisotropy 



[2 0l | as well as the helicity distribution typical for pseu- 
doscalar Dalitz decays, although attenuated. Figure [T4l 
shows the ratio of the reconstructed/simulated center-of- 
mass polar {dN / dO c . m ) and helicity (dN/da) distribu- 
tions for three different pair mass bins. Note that these 
ratios have been reflected about 90° and both halves 
added in order to reduce statistical fluctuations. The 
normalization is arbitrary. The resulting angular dis- 
tributions exhibit anisotropies which are quantified by 
adjusting 1 + A 2 cos 2 0c. m. and 1 ± B cos^ a forms, re- 
spectively. 

The low-mass anisotropies are large and consistent 
with our expectations for the neutral pion. The fitted 
polar coefficient A 2 — 0.61 ±0.09 corresponds, according 
to our simulation, to an un-attenuated A 2 = 0.76 ±0.11, 
in agreement with the polar anisotropies of charged pions 
observed in Ar+KCl (20|, namely A 2 = 0.75 ± 0.05. The 
helicity, B — 0.71±0.05, is attenuated by the thermal 
emission of the pion from its QED value B = 1, again 
consistent with the expectation from our simulations. 

Intermediate-mass pairs are more interesting because 
only about 25 - 30 % of their yield is exhausted by -q 
Dalitz pairs, the dominant excess part being of non- 
trivial nature (see discussion in III.B). Angular distri- 
butions might provide some constraints on its possible 
composition. The intermediate-mass bin in Fig. [14] dis- 
plays large anisotropies as well, both for polar emission 
angles, with A 2 = 0.72 ± 0.24, and for helicity, with 
B = 0.55 ± 0.12. Taking into account the attenua- 
tion, this is again compatible with the 1 + cos 2 a be- 
havior typical for pseudoscalar meson, but also A de- 
cays. Subtracting the simulated rj contribution from the 
data, the pure excess angular distributions have been 
obtained. They are represented as well in the figure, 
together with the corresponding fits, showing that the 
anisotropies (Af = 0.69 ± 0.30 and B x = 0.51 ± 0.17) of 
the excess yield turn out to be very similar to the ones 
of the rj. This suggests that a large fraction of the ex- 
cess can be attributed to decays of the A resonance for 
which indeed B — 1 is also expected. One has to keep in 
mind, however, that the nucleon-nucleon bremsstrahlung 
contributes as well in this mass region and has to be con- 
sidered in a full description. 

The high-mass bin is unfortunately very low in statis- 
tics, but exhibits nevertheless a strong polar anisotropy 
(A 2 = 2.2±1.5), whereas its helicity distribution is within 
its (large) statistical errors compatible with B = 0. This 
is to be expected, if this mass bin contains large contri- 
butions from vector mesons emitted from a completely 
thermalized source, just like it was observed in the NA60 
experiment [53j |. although this is in conflict with the ob- 
served large slope parameters. 



IV. CONCLUSIONS AND OUTLOOK 



The results on e + e production obtained with HADES 
in the medium- heavy 40 Ar+KCl system at 1.756A GeV 
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FIG. 14: (Color online) Ratio of measured and simulated dielectron center-of-mass polar distributions dN /d0 c .m. (a) and helicity 
distributions dN/da (b) for three mass bins: M ee < 0.15, 0.15 < M ee < 0.5, and Af ee > 0.5 GeV/c 2 (from top to bottom). 
The error bars are statistical. The Pluto cocktail simulation was done assuming isotropic emission and decay of the dileptons 
(see text). The curves are fits to the data yielding the anisotropy coefficients, i.e. polar A2 and helicity B. The coefficients A% 
and B x result from fits (dashed curves) to the ^-subtracted ratios (triangles). 



show an intermediate-mass pair excess over long-lived 
sources a factor 2-3 stronger than the one observed in 
elementary and 12 C+ 12 C reactions. We have discussed 
the enhancement in the integral pair yields as well as 



differentially in the pair mass and transverse mass distri- 
butions. We have argued that this behavior signals the 
onset of influence of the nuclear medium on dilepton pro- 
duction. We have isolated the emission from the medium 
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by subtracting the known contributions from long-lived 
pair sources. By presenting transverse mass and angular 
distributions of the eta-subtracted yields, we have been 
able to characterize this contribution further. From the 
analysis of the excess transverse-mass slope and angu- 
lar anisotropies we concluded that they are compatible 
largely with A Dalitz decays, suggestive of resonance 
matter. 

Furthermore, for the first time at SIS18 energies, a 
clear ui signal could be observed in heavy-ion collisions, 
quantified and compared with the prediction of a statis- 
tical hadronization model. This result allows, in particu- 
lar, to put tight constraints on vector meson production 
in heavy-ion collisions at beam energies of a few GeV. 
From the shape and position of the observed peak no 
direct indications could, however, be found for medium 
modifications of this vector meson. In fact, in case of 
the very strong broadening of the ui implied by the in- 
terpretation of photoproduction data, our measurement 
would anyhow have been sensitive only to the freeze-out, 
i.e. vacuum part of its yield. 

A first and preliminary comparison of our Ar+KCl in- 
variant mass e + e~ spectrum with predictions of models 
based on transport theory has been discussed in [ill ■ 
While, both the Hadron String Dynamics (HSD) model 
[72| and the Ultra-relativistic Quantum Molecular Dy- 
namics (UrQMD) model [73[ achieve good agreement for 
invariant masses below 0.15 GeV/c 2 , at intermediate and 
higher masses the description of the pair yield is not 
yet satisfactory. Known reasons are the still imperfect 
description of some of the elementary processes imple- 
mented in transport models, namely bremsstrahlung and 
vector meson production [18]], but also the largely open 
question about how to treat possible in-medium modi- 
fications of these processes. On the other hand, these 
models suggest that the part of our spectrum most sen- 
sitive to possible in-medium modifications should be the 
region of excess yield, namely the dielectrons with masses 
of 0.5-0.8 GeV/c 2 , which hence need to be characterized 
in detail. With the additional and more differential e + e~ 
data presented in this paper we have provided the infor- 
mation required to improve the theoretical description 
of dilepton production in heavy-ion collisions. Further- 
more, once supplemented with data on yet heavier re- 
action systems, these new results are expected to reveal 
and quantify medium modifications of hadrons in warm 
and dense nuclear matter. 

In summary, we have presented phase space distribu- 
tions (invariant mass, transverse mass, rapidity, polar 
angle, and helicity angle) of dielectrons for the reaction 



Ar+KCl at 1.756A GeV. From these the following results 
have been obtained: (i) observation of a strong excess 
(up to factor 3) over N+N collisions of the pair yield at 
intermediate masses; (ii) first observation of uj produc- 
tion in heavy-ion collisions at such a low beam energy, 
yielding a <j)/oj ratio consistent with maximal violation of 
OZI suppression; (iii) a systematic study of transverse- 
mass distributions with the observation of unexpectedly 
large inverse-slope parameters of up to 135 MeV, which 
might be related to final-state processes and/or the spec- 
tral functions of the contributing dilepton sources; and 
(iv) first exploitation of the virtual photon polarization 
observable at low beam energies. 

Our studies on dielectron production in heavy-ion re- 
actions will be pursued over the next years with an up- 
graded HADES detector which will have the ability 
to handle the high track densities from truly heavy colli- 
sions, in particular also the 197 Au+ 197 Au system. These 
data runs will thus provide the full systematics required 
to address open questions on the origin and properties 
of the intermediate-mass pair yield. In parallel to this 
heavy-ion program, the HADES experiment will also take 
up studies making use of the GSI secondary pion beams, 
in elementary 7r+p reactions as well as in 7r + A reactions. 
The pion-beam experiments will allow to conduct a com- 
prehensive study of the contribution of specific baryon 
resonances to dielectron emission, in vacuum and in the 
nuclear medium. The information gained this way will 
in turn also help to improve our understanding of dilep- 
ton radiation from the hot and dense hadronic medium 
produced in the reactions with heavy ions. 
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